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Abstract 

The recently proposed spin-aligned neutron-proton pair coupling scheme is 
studied within a non-orthogonal basis in term of the multistep shell model. 
This allows us to identify simultaneously the roles played by other configu- 
rations such as the normal pairing term. The model is applied to four-, six- 
and eight-hole N = Z nuclei below the core ^°°Sn. 
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Many features in nuclear structure physics can be understood in term of 
the seniority coupling scheme, which was first introduced in atomic physics 
by Racah [l| . This scheme showed to be extremely useful for the classification 
of nuclear states in the jj-scheme , particularly in semimagic nuclei 

with only one type of nucleons. The lowest-seniority pair (with v = 0) has 
nothing special from a coupling point of view since the nuclear state can 
then be constructed in a variety of equivalent ways through other pairs. In 
particular, the aligned like-nucleon pair coupling was proposed in Ref. 5|, 
which may manifest itself from the energy differences of mirror nuclei 6|. 
The driving force behind the dominance of seniority coupling is the strong 
pairing interaction between like particles. 

The neutron-proton (np) correlation breaks the seniority symmetry in 
a major way. Correspondingly, the wave function is a mixture of many 
components with different seniority quantum numbers. It is not clear yet how 
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this kind of states can be classified in the jj-scheme. The stretch scheme, 
which corresponds to the maximally aligned intrinsic angular momentum, 
was proposed in the 1960s to describe the rotational-like spectra of open- 
shell nuclei 0, Isj . But now it is widely accepted that a proper description of 
deformation involves the mixture of different orbitals. 

The low-lying yrast states in IgPd were recently reported [9| . This is the 
heaviest N = Z nucleus with measured spectrum so far. It was suggested 
that in this nucleus, as well as in neighboring nuclei like ^^Cd, the properties 



of the low-lying states can be largely described in a single Ogg/2 shell |10l. Ill 
Furthermore, it was proposed that the low-lying yrast states in these N = Z 
nuclei can be classified by a spin-aligned np pair coupling scheme jo], [lo| . 
That is, the ground state wave functions do not consist mainly of pairs of 
neutrons [uu) and protons (tttt) coupled to zero angular momenta, but rather 
of isoscalar np pairs (z/vr) coupled to the maximum angular momentum J, 
which in the shell OgQ/2 is J = 9 10|. A detailed shell-model analysis 



of the spin-aligned np pair coupling was performed in Ref. [10[. This shell 
model calculation was done in the standard fashion of using as representation 
the tensorial product of neutron times proton degrees of freedom, thus easily 
takes into account the Pauli principle. The drawback with this representation 
is that it is not straightforward to realize that the states are mainly deter- 
mined by np pair degrees of freedom. This can only be done by projecting 
the shell model wave function into the particular np component one wishes 
by using two-particle coefficients of fractional parentage. This calculation 
becomes rather involved for systems with more than three pairs. Moreover, 
it does not allow one to study simultaneously the competing effects of differ- 
ent np pairs, which requires a large amount of independent projections. A 



similar calculation was done in Ref. [ll| for the two-pair case, just confirm- 



ing the coupling scheme mentioned above. For systems with three and four 
pairs, the interacting boson model was applied through the boson mapping 



of the aligned np pair [11 



In this Letter we will show that a suitable representation to study simul- 
taneously different partitions of a system consisting of many np excitations 
is the multistep shell model method (MSM) [l2|. In this method one solves 
the shell-model equation in several steps. In the first step one constructs the 
two-particle states. In the second step one proceed by solving the three- or 
four-particle states in terms of the two-particle states calculated in the first 
step. In our case we will solve the two-neutron plus two-proton system within 
a non-orthogonal overcomplete basis in terms of the (z/vr) (z^tt) excitations 
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at the same time as the (uu) (tttt) ones. With the four-particle system thus 
evaluated, we will proceed to evaluate the six-particle system in terms of the 
coupling of the four-particle and two-particle states. For the eight-particle 
system one can choose the MSM basis such that it consists of the products 
of the four-particle states in the form (z/tt) (^^tt). Systems with more pairs 
can be described in the same fashion in successive steps. 

The MSM treatment of four and six identical particles was performed in 
Ref. [12)]. Although the formalism to be used here is similar, the present 
calculation is even more challenging due to the presence of both neutron 
and proton degrees of freedom. In the case of identical particles the MSM 
basis is overcomplete mainly because it violate the Pauli principle. In our 
case the overcompleteness of the basis is even more severe since our basis 
elements may count twice the same states besides violations of the Pauli 
principle. For instance, in the two-pair case the basis elements {i>i/)^(jrn) and 



(i^tt) ® (z^vr) may be proportional to each other (see, also, Refs. [10|, llli)- The 



overcounting thus occurring is a result of describing the np and like-particle 
excitations at the same time. In the MSM this complication is overcome by 
evaluating the overlap matrix from which an orthonormal set of states can 
be constructed. This can be a very time consuming procedure for systems 
with more than two pairs. 

We will use the Greek letter 7„ to label the n-particle np states. Since 
we will only consider cases with equal number of neutrons and protons out- 
side a closed shell, n will be an even number. A m-proton (m-neutron) 
state will be labelled by am (/9m) • Therefore the np states will be I72) = 
-P"'' (72)10) where the np creation operator is ^"'"(72) = J2ip-^iw'^l2)clcp 
and cf (Cp ) is the neutron (proton) single-particle creation operator. In 
the same fashion the two-proton (two-neutron) creation operator will be de- 
noted as P+(a2) (P+(/32)) (c.f., Eq. (9) of Ref. 0). The four-particle state, 
|74)=P+(74)|0), is 

P+(74) = J]X(«2/32;74)P+(«2)P+(/32) 

02,^2 

+ 5^ X(7272;74)P+(72)P+(72), (1) 

72<72 

where all possible like-particle and np pairs are taken into account. Since the 
number of MSM basis vectors is larger than the dimension of the shell model 
space, the wave function amplitudes X are not well defined in our case and. 
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therefore, they are not meaningful physically. The meaningful quantities are 
the projections of the basis vectors upon the physical vector, which we denote 
as 



F(a2/32;74) = (74|P+(a2)P+(/32)|0), 

Fil2l2-n4) = (74|P+(72)P+(72)|0). (2) 

The orthonormality condition now reads 

Kii = ^("2/32; 74)^^(02/32; 74) 

0:2,132 

+ ^ X(727^;74)F(727^;74). (3) 

72<72 

The norm of the MSM basis 1727^) = P+(72)P+(7^)|0), i.e., Ar(7272;74) = 
V (727217272)5 iiot be unity. Therefore the interesting quantity is not 
the projection F but rather the cosine of the angle between the basis vector 
and the physical vector, i.e., cos(0) = x and 

2;(7272; 74) = ^(7272; 74)/A^(7272; 74)- (4) 

If we would have taken as basis elements the complete set of orthonormal 
states {P~^a2)P~^{(32)\0)} (which is the standard shell model basis as used 
in Ref. |10|) then the second term in Eq. ([3]) would not have appeared 
and one would have obtained X {a2(32] 'Ja) = a;*(a2/32; 74), as expected in an 
orthonormal basis. One thus sees that the advantage of the MSM basis is 
that one can extract the physical structure of the calculated states just by 
examining the quantity x. 

For the six-particle case we will use the MSM partition of two- times 



four-particles, as it was done in Ref. [12| for systems with six like particles. 



Thus the corresponding wave function will be I76) = P"'"(76)|0), where 

^^(76) = Yl ^(^274; 76)P+(72)P+(74)- (5) 



72,74 



As before, we will evaluate the projection of the basis vectors upon the phys- 
ical vectors, i.e., -^(7274; 76); and the corresponding cosine function x. In 
this six-particle case one can also view the MSM basis elements as the direct 
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tensorial product of three pairs. This is a unique feature of the MSM. The 
projection of such a MSM basis upon the physical vector is, 



i^(727^72'; 76) = (76|P+(72)P+(7^)P+(7D|0) 

= X^^(7274;76)i^(7272;74), (6) 

74 

from which one can evaluate the norm as A^^(727272 j ^) = X]76(a) 1-^(727272 j 76) P 
where A is the total angular momentum of the state and the sum runs over 
all physical states 76 with angular momentum A. The cosine of the angle 
between a MSM basis and the physical state is 

a:(727272 ; 76) = i"(727272 ; 76) /iV (727272 ; A). (7) 

We will describe the eight-particle states as |78) = P~'"(78)|0), where 
-P^(78) = Z]74<74-^(7474;78)-P"^(74)-P"^(74)- Proceeding as above we will 
also evaluate the cosine of the angle between |78) and all the possible four- 
pair states that can be formed. 

It is important to point out that for any MSM basis element corre- 
sponding to the n-particle system it is a;^(6„, 7„)=1. This is because the 
vectors |7„), which are eigenvectors of the n-particle Shell Model Hamilto- 
nian, form an orthonormal (complete) set. Therefore the cosine x(6„,7„) is 
the probability of the state |7„) occupying the basis state |6„). 

We will apply the method to study the spin-aligned np pair coupling 
scheme j^, [lO]. We will restrict our calculations to the single Ogg/2 shell 
with the interaction matrix elements taken from Ref. [lo|. But it should be 
emphasized that the formalism proposed in the present work can be naturally 
generalized to systems with many shells. 

In the cases of ^^Cd and ^^Pd the low-lying spectra are determined by 
the isoscalar and strongly attractive matrix element ((fi'9/2)^; 9|V^|((79/2)^; 9) 
Toll, which corresponds to the maximally aligned np pair configuration. The 



extend to which this determines the spectrum can be deemed by the evolution 
of the calculated levels as a function of the variations of that matrix element. 
Calling Vg{6) = Vg{0){l + 6) it is found that as Vq{5) — )• the spectrum 
tends to have a seniority-like form. At the other extreme, approaching 5=1, 



a tendency towards a vibrational-like spectrum seems to take place [10| . This 
is not surprising since, as also seen below, all low-lying yrast states in both 
nuclei are isoscalar np pair excitations |10| . 
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Figure 1: (Color online) Upper: Coefficients corresponding to the wave function of 
the ground state of ^^Cd as a function of the controlling parameter 5 (see text). The 
labels indicate 1: \[[vTT)gf]Q); 2: |(zyi/)o(7r7r)o; 0); 3: |((i/7r)i)2; 0); 4: |(z/i/)2(7r7r)2; 0). 
Lower: Same as the upper panel but for the first 2+ state. The labels correspond to 
1: |((i^7r)9)^; 2); 2: \{vtt)q{vti)2:2) . Only configurations with > 0.4 are shown for 
simplicity. 

The MSM provides in a straightforward fashion the structure of the states 
in terms of all possible configurations. Thus, in Fig. [1] we show the main val- 
ues of the probabilities (Eq. (jl])), as a function of the controlling parameter 
6, for the ground state and first 2+ state of ^^Cd. The striking feature in 
this figure is that the spectrum for 5=0 is dominated by the isoscalar config- 
uration ((z/7r)9)^. Moreover, one sees that as Vg{6) becomes more attractive 
(5 — )• 1) the pairing state becomes less and less relevant while the importance 
of other isoscalar components increases. This feature is even more remarkable 
for the state 2f, where the spin-aligned np pair coupling ((z/7r)9)^ dominates 
the wave function for all values of S shown in the figure. 

One may wonder whether the pairing mode |(z/i/)o(7r7r)o; 0), which does 
not dominate the ground state in ^^Cd, would be located as a state 0+ higher 
up in the spectrum. However, this is not the case, as can be seen from Tabled] 
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Table 1: Leading configurations in the first five states of ^^Cd for a given total angular 
momentum /. 





7 = 




7 = 2 




7 = 4 




7 = 6 




n 


Configuration 


\x\ 


Configuration 


|.x| 


Configiuation 


|x| 


Configuration 


\x\ 


1 




0.96 


(l^7r)g (S) (!^7r)9 


0.99 




0.97 




0.84 


2 


{I'n)^ (S) (i^7r)5 


0.88 


(i/7r)8 (gi (i^7r)9 


0.92 


{u7t)s (g> (1^77)9 


0.93 


(1/71)7 CS) {i/7r)g 


0.60 


3 


(v7r)o (S) (i^7r)o 


0.77 


(i/7r)7 (g) (/^7r)9 


0.94 


(j'7r)7 (S (/^7r)9 


0.90 


(i'!^)o(6) "S (7r7r)6(o) 


0.94 


4 


(S) (i^Tr):; 


0.91 


(i/7r)7 (g {un)s 


0.82 


(i/7r)5 (g> (;^7r)9 


0.90 


(1^1^)2(6) "S (7r7r)6(2) 


0.65 


5 


(l'7r)4 (S) 


0.80 


{utt)^ (Si (;^7r)7 


0.64 


(i/7r)6 IS (/^7r)9 


0.92 


(i^f )2{6) "S ('r7r)6(2) 


0.84 



where we listed the main components of the lowest-lying states in ^^Cd for dif- 
ferent total angular momenta. Rather it is distributed throughout the spec- 
trum. In contrast, the isoscalar aligned mode |((z/7r)9)^; /) is mainly concen- 
trated in the yrast states. It may seem weird that these two modes produce 
different results since they may be related to each other just by an exchange 
of neutrons and protons. But they are not the same, as shown by the angle (j) 
between them which gives cos(0)=O.62, i.e., = 52°. It is also interesting to 
notice that the norm of the aligned state is (((z/7r)9)^; 0| ((1/71)9)^; 0) = 2.00001, 
which shows that the influence of the Pauli principle upon |((z/7r)9)^; 0) is neg- 
ligible and, therefore, it represents virtually a bosonic mode (see also Ref. 

H). 

The calculated spectrum of the odd-odd nucleus ^^Ag is shown in Fig. |2l 
where we have grouped the levels according to their isospin. There are only 



two states which have been measured in this case [13|, namely the O"*" ground 
state and a state 21+ at 6.67 MeV. These states are in reasonable agreement 
with experiment and coincide with previous shell model calculations ij]. The 
strong influence of the aligned isoscalar matrix element Vg upon the T = 
states can be inferred from the figure, where it is seen that the energies of 
T = states are much more sensitive to the controlling parameter 6 than 
those with T = 1. 

The T = states in ^^Ag are specially interesting because in this case the 
MSM basis vector consisting of the three 9"*" isoscalar aligned states is not 
hindered by any symmetry (recall that only states with total isospin T = 
can be coupled from the isoscalar np pairs). Indeed we found that most of 
the yrast levels in Fig. [2l except the I61 and I81 states, are mainly built 
by the isoscalar aligned np pairs. As an example, in the lower panel of Fig. 
I2]we present the the probabilities for main components of the 7f state of 
^"^Ag, which is calculated to be the lowest T = state, as a function of the 
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Figure 2: (Color online) Shell model spectra of ^^Ag calculated in the 6(79/2 shell as a 
function of the controlling parameter S. 
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Figure 3: (Color online) Upper: Coefficients corresponding to the wave func- 
tion of 94 Ag (Of) as a function of 6. The labels indicate 1: |(z/7r)o(i^7r)9(i^7r)9); 0); 
2: |(zy7r)o(H2(^^)2);0); 3: \{v^\(vv)^{-ki,\)-Q)- 4: |(i.^)o(i^7r)o(z^7r)o); 0); 5: 
(z/7r)o(i^7r)5(i^7r)5); 0); 6: |(i'7r)o(j^i')o(7''7r)o); 0). Lower: Same as the upper panel but 
for the l\ state. The labels indicate 1: |(zy7r)o(i^7r)9(z^7r)9); 7); 2: |(z^7r)7(zyi/)o(7r7r)o); 0); 3: 
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controlling parameter 5. 

It is seen from Fig. |2]that at 5 = the ground state of ^^Ag carries T = 1. 
It may thus seem that in odd-odd system the isovector pairing mode retakes 
its predominance. We found that this is not the case by analyzing the ground 
state wave function in terms of the tensor product of three pairs, as can be 
seen from the upper panel of Fig. [31 The most important configuration con- 
sists of the spin-aligned np pair state |(i/7r)9(i/7r)9); 0) that was dominant in 
^^Cd and the other possible configuration, i.e., the isovector state |(z/7r)o). In 
fact the low-lying T = 1 states here have the same origin as the correspond- 
ing T = 1 states in ^^Cd. One sees that in this odd-odd six-particle case the 
probabilities of the ground state occupying different basis states are larger 
than in ^^Cd, Fig. [T] This is not surprising since for systems with more than 
two pairs there are many nearly equivalent combinations that can be built 
in the same fashion. 

In the analysis of eight-particle systems like ^^Pd we choose as MSM basis 
the partition |7474;78)- In this case the MSM basis is highly overcomplete. 
For instance there are 36 shel-model O"*" states while the corresponding MSM 
dimension is 915. Within this basis we calculated in Table |5] the quantities 
X, i.e., the cosines of the angles between the vectors I78) and all the possi- 
ble vectors that can be formed by the coupling of four pairs. Since many 
combinations are similar to each other there is not a value of x which is sig- 
nificantly larger than the others. But one finds, again, that for the ground 
state of ^^Pd the most important MSM configuration is the one correspond- 
ing to the four aligned pairs. The second one is a combination of two 
aligned 9+ states and the normal pairing states. This is expected since in 
the two-pair case of ^^Cd the second largest component is the normal pairing 
term. An important feature in this case is that for the pairing state it is 
x^{a2 = 0+,/32 = O+a'2 = 0+/3^ = 0+;78 = 0+) = 0.46. This is a relatively 
small number. Indeed it occupies the 10th place in order of importance. 
This reflects, once again, the dominance of the aligned configuration in this 
nuclear region. 

Summarizing, we have in this paper extended the MSM method proposed 



in Ref. [12[ to incorporate both neutron and proton degrees of freedom 
and applied it to study the recently proposed spin-aligned np pair coupling 
scheme. We have applied the method to analyze four-, six- and eight-hole 
states in N = Z nuclei below the core ^°°Sn. The calculations were performed 
within the restricted Ogg/2 shell for simplicity. But this work opens the way 
for even more challenging calculations, involving many particles and/or many 
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Tabic 2: Configurations witli tlie largest probabilities for the state ^^Pd(0]^) corresponding 
to the tensorial products of different two-particle states (upper) and four-particle states 
(lower). 



Configuration 


x'^ 


\l2 


= 9+7^ = 9+7^' = 9+7^" 


= 9+) 


0.85 


\l2 


= 9+7^ = 9+^2 = 0+/32 


= 0+) 


0.76 


72 


= 8+7^ = l+a2 = 0+^2 


= 8+) 


0.56 


\l2 


= 8+7^ = l+a2 = 8+/52 


= 0+) 


0.56 


\l2 


= 1+7^ = l+a2 = 0+/32 


= 0+) 


0.52 



shells. It would also allow one to truncate the shell model basis in terms of 
spin-aligned np pairs or other coupling schemes. 
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